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Abstract H l observations of high-velocity clouds (HVCs) in- 
dicate, that they are interacting with their ambient medium. 
Even clouds located in the very outer Galactic halo or the in- 
tergalactic space seem to interact with their ambient medium. 
In this paper, we investigate the dynamical evolution of high 
velocity neutral gas clouds moving through a hot magnetized 
ambient plasma by means of two-dimensional magnetohydro- 
dynamic plasma-neutral gas simulations. This situation is rep- 
resentative for the fast moving dense neutral gas cloudlets in 
the Magellanic Stream as well as for high velocity clouds in 
^S) ' general. The question on the dynamical and thermal stabiliza- 
' tion of a cold dense neutral cloud in a hot thin ambient halo 
plasma is numerically investigated. The simulations show the 
formation of a comet-like head-tail structure combined with a 
, magnetic barrier of increased field strength which exerts a sta- 
1-^ bilizing pressure on the cloud and hinders hot plasma from dif- 
Qv^ fusing into the cloud. The simulations can explain both the 
I ■ survival times in the intergalactic medium and the existence of 
, head- tail high velocity clouds. 

Key words: Magnetohydrodynamics (MHD) - Galaxies: ha- 
los - intergalactic medium - Magellanic Clouds - Galaxies: 
magnetic fields 
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1. Introduction 

High-velocity clouds (HVCs) - first discovered by MuUer et al. 
( 1963 ) - are defined as neutral atomic hydrogen clouds with ra- 



dial velocities that cannot be explained by simple galactic rota- 
tion models. After almost 40 years of eager investigations there 
is still no general consensus on the origin and the basic physical 
parameters of HVCs (e.g. Bregman 1999 ). The most critical is- 
sueof HVC research i s the d istance uncertainty . Dan ly et al. 
( |1993| ), Keenan et al. ( |1995| ) and Ryans et al. ( |1997| ) consis- 
tently determined an upper distance limit of d < 5 kpc to HVC 
complex M. The most important step forward is the very re- 
cently determined distance bracket of 4 < d < 10 kpc towards 
HVC complex A by van Woerden et al. (1999). These results 



clearly place the HVC complexes M and A in the gaseous halo 
of the Milky Way. 

Parallel to the growing evidence that a significant fraction 
of the HVC c omplexes are located in the Milky Way halo. Blitz 
et al. (1999) supported the hypothesis that some HVCs are 
of extragalactic origin. They argued, that it is reasonable to 
assume that primordial gas - left over from the formation of 
the Loc al Gro up galaxies - may appear as HVCs. Braun & 
Burton (1999) identified 65 compact and isolated HVCs and 
argued that this ensemble represents a homogeneous subsam- 
ple of HVCs at extragalactic distances. Observational evi- 
dence for extragalactic HVCs may also be found by the de- 
tection of th e high ly ionized high-velocity gas clouds by Sem- 
bach et al. ( |1999| ), because of its very low pressure of about p 
k-i « 5 Kcm^ 

The Magellanic Stream (MS) and the Leading Arm (LA) 
(Putman et al. 1998) both form coherent structures over sev- 



eral tens of degrees having radial velocities in the HVC regime. 
They represent debris most likely caused by the tidal interac- 
tion of the Magellanic Clouds with the Milky Way. The dis- 
tances to these features are of the order of 50 kpc. 



Meyerdierks (1991) detected a HVC that appears like a 
cometary shaped cloud with a central core and an asymmet- 
ric envelope of warm neutral atomic hydrogen (the particular 



HVC is denoted in literature as HVC A2). He interpreted this 
head-tail structure as the result of an interaction between the 
HVC and normal galactic gas at lower velocities. Towards the 
HVC complex C, Pietz et al. (1996) discovered the so-called 



Hi "velocity bridges" which seem to connect the HVCs with 
the normal rotating interstellar medium. The most straight for- 
ward interpretation for the existence of such structures is to 
assume that a fraction of the HV C gas was stripped off the 
main condensation. Briin s et a l. (200C) extended the inves- 
tigations of Meyerdierks (1991) and Pietz et al. (1996) over 
the entire sky covered by the new Le iden/D wingeloo Hi 21- 
cm line survey (Hartmann & Burton 1997) and found head- 
tail structures in all HVC-complexes including the Magellanic 
Stream, except for the very faint HVC-complex L. Their anal- 
ysis revealed that the absolute value of the radial velocity of 
the tail is always lower than the value for the head of the HVC 
(|wLSR,taii| < IwLSR.hoadD- In addition, it was shown that the 
fraction of HVCs showing a head-tail structure increases pro- 
portional to the peak column density and increasing radial ve- 
locity |wgsr|- 

HVCs mostly appear as "pure" neutral atomic hydrogen 
clouds. Absorption line studies provide information on the ion- 
ization state and the metalicity of HVCs. The results indicate 
that the bright and very extended HVC complexes consist (at 
least partly) of processed material, having < 1/3 of the solar 
abundances (Wakker 2001 ). 

Thus, the question of the dynamical and thermal stabilization of 
a cold dense neutral cloud in a hot thin ambient plasma arises. 
To confine a HVC in the Magellanic Stream by pressure arel 
atively high halo gas density is required (Mirabel et al. 1979 
Weiner & Williams |1996| ). Therefore, the lifetime of the HVCs 
should be significantly limited by evaporation (Murali 200C| ) 



since the necessary pressure for confinement is associated with 
a high energy transfer. An alternative mechanism of dynami- 
cal stabilization is the magnetic confinement of the cloud. In 
this contribution, we numerically investigate the formation of 
a magnetic barrier around the HVC in the Magellanic Stream 
and its stabilizing effects on the neutral gas cloud. 
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Figure 1. This figure shows the column density distribution of 
a part of the Magellanic Stream observed with the multi-beam 
facility of the Parkes telescope. The contour levels represent 
0.1, 1,2, 3, ...-lOi^cm-^. 



2. Observations 

The gaseous arms of the Magellanic System are very inter- 
esting objects to study the evolution of clouds moving with 
high velocities relative to an ambient medium. The Magel- 
lanic Stream (MS) is a coherent structure which starts from 
the Magellanic Clouds and trails over 100° passing the south- 
ern Galactic pole. The physical parameters of both the gas in 
the Magellanic Stream and the ambient medium must support 
survival times of the order 10^ yr predicted by numerical sim- 



ulations (Gardiner 1999 1 



The low distance of about 50 kpc enables studies of the 
distribution of gas on scales of about 200 pc with single dish 
telescopes. Recently, the Parkes multi-beam narrow-band sur- 
vey of the entire Magellanic System was completed (see Briins 
et al. 2000). Fig. |l| shows exemplary the column density distri- 
bution of a part of the MS, also known as MS IV. The Stream is 
subdivided in two filaments containing a number of condensa- 
tions. Some of them show a horseshoe like structure. In addi- 
tion, the Stream is accompanied by a number of isolated clouds 
at similar velocities that are most likely part of the gaseous fea- 
ture. The detailed analysis of the gas within the MS is difficult, 
as some clouds are superposed on the same line of sight. The 
analysis of the isolated clouds is much easier in this respect. 
Observable quantities are the angular extent, the line-widths, 
and the column density distribution. Other important values 
like the gaseous mass, the density etc. can only be calculated. 



if the distance to the clouds is known. We assume a distance of 
50 kpc to all clouds in this section. 

Typical line-widths derived from isolated clouds are 
FWHM « 25 - 30 km s^^, typical diameters are of the or- 
der 40'. This corresponds to 580 pc at a distance of 50 kpc. 
A particle moving with a velocity of v w 15 km s^^ w 15 
kpc Gyr^^ crosses the cloud in « 0.04 Gyr, i.e. the crossing 
time is much smaller than the age of the Magellanic Stream. 
The virial mass of a cloud can be calculated using A/vii = 
210 r (Av)^, where r is the radius of the cloud in pc and Av is 
the linewidth in km s^^. Virial masses for clouds in the MS are 
Mvir = 3.94 10^ (r / 300pc) (Av / 25km ^-^f Mq, while ob- 
served HI masses are of the order few x 10^ Mq. Apparently, 
the observed HI masses are orders of magnitude to low to pro- 
vide gravitational stabilization. Consequently, the cloud must 
be confined by external pressure. Direct observations of the 
ambient medium provide only weak constraints on the physi- 
cal parameters. The density of the ambient medium is of the 



order tt-igm ~ 10 ^ cm (Murali 2000; Weiner & Williams 



1996; Mirabel et al. 1979). The temperature of the gas in the 



outer Galactic halo is expected to be similar to the tempera- 
ture in the lower Galactic halo (Tigm ~ 1 - 2 x 10^ K). For 
a fully ionized medium, the pressure can be calculated using 
PiGM k"^ = 2.3 niGM Tigm, i-e. Pigm k"^ w 50 K cm"^. 

The pressure within the MS can be estimated using the cor- 
relation between linewidth, turbulent velocity and kinetic Tem- 
perature 



Aw = 2 ln2 



f2kTk 



1/2 



where Av is the linewidth, Tk is the kinetic temperature, k is 
the Boltzmann constant, mn is the mass of the hydrogen atom, 
and V is the most probable turbulent velocity. The assumption 
that there is no turbulence provides a possibility to derive an 
upper limit to the kinetic temperature (a.k.a. Doppler temper- 
ature) from Eq. To = 21.8 (Avf « 2x10* K, where 
Td is in K and Av is FWHM in km s"^ The actual ki- 
netic temperature must be lower, as processes like turbulence 
are likely to occur Typical observed volume densities within 
the Magellanic Stream are of the order n = few x I0~^ cm^'', 
assuming a distance of 50 kpc for the Magellanic Stream. Us- 
ing Pk^^ = nT allows to estimate the pressure in the Stream 
to be of the order Pk~^ ~ few x 100 Kcm""^. 

Pressure equilibrium between the neutral gas in the clouds 
and the ionized medium is possible, as both pressures are only 
rough estimates. In addition, the contribution of magnetic 
fields to the question of stability is completely unknown. 

There are a number of clouds in the Magellanic Stream 
and the Leading Arm showing up with head-tail structures, i.e. 
asymmetries in the column density distribution associated with 
velocity gradients (Briins et al. 200Cj ). A detailed analysis of 
these features is extremely complicated, because of the super- 
position of clouds along the same line of sight. 

Another interesting class of HVCs are the so-called com- 



(1) 



pact, isolated HVCs. Braun & Burton (1999) argue, that these 



HVCs are located within the intergalactic medium of the Lo- 
cal Group. They are by definition compact and isolated - a 
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Figure 2. This figure shows the column density distribution of 
HVCl 25+4 1-207. Contour Hnes represent 1, 2, 4, 10, 20, 40, 
100 X 10^* cm^^. This clouds shows a very prominent head- 
tail structure. 



detailed analysis of the HI distribution is therefore much eas- 
ier compared with the Magellanic Stream. Fig. ^ shows the 
column density distribution of a compact and isolated high- 
velocity cloud with a distinct head-tail structure. This cloud 
was discussed in detail in Briins et al. (2001 1. The distance to 
this cloud was determined to be of the order 130 kpc, i.e. the 
cloud is located in the very outer halo of the Milky Way or the 
near intergalactic space. We consider this cloud as a generic 
example for the class of head-tail HVCs investigated in this 
contribution. 

Fig. H shows the result from a Gaussian decomposition of 
the HI data of HVC125+41-207. About 12.5 % of the Hl-gas 
is located in the tail. The spectra show in general two compo- 
nents, one with a low and one with a high velocity dispersion 
(i.e. a cold and a warm gas-phase). The velocity dispersion of 
the warm Hl-gas is increasing from FWHM w 10 km s^^ to 
FWHM w 20 km s^^ over the extend of the cloud and re- 
mains approximately constant in the tail. This corresponds to 
a Doppler temperature increasing from Td = 2000 K to Td = 
10000 K in the tail. There is also a trend of higher velocity 
dispersions towards the outskirts of the cloud. In addition, the 
radial velocities |wgsr| of the warm component are lower than 
the velocities in the colder gas. This is a strong indication that 
the warm gas is in the process of being stripped off the main 
body of the HVC. The observed morphology of HVC125+41- 
207 clearly demonstrates, that this cloud is currently interacting 
with its ambient medium similar to the HVCs in the Magellanic 
Stream and Leading Arm. 

Briins (in prep.) performed deep HI observations of a sam- 
ple of 10 compact HVCs with the 100-m Effelsberg telescope. 
All of them show up with a head-tail structure - some of the 
tails have a fairly low column density (wlO^® cm~^). 

Altogether, the observations demonstrate, that HVCs in the 
Galactic halo or in the near intergalactic space interact with 
their ambient medium producing head-tail structures. 
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Figure 3. This figure shows the result from a Gaussian decom- 



position of the HI data of HVC 125+4 1-207 (Bruns et al. |2001[ ). 
The cold gas-phase is represented by filled circles, the warm 
gas-phase by open circles. Results from a single Gaussian fit 
are represented by crosses. The diagrams a,c and e show the 
results for a slice along the tail, the diagrams b,d and f repre- 
sent the slice perpendicular to the orientation of the tail. The 
line in diagram a is an exponential fit to the data. The line in 
diagram e is a fit to the velocity field. 



3. Numerical simulations of high velocity clouds 

In order to study the dynamical evolution of a fast moving, 
cool, dense neutral gas cloud in a magnetized, hot ambient 
medium we performed 2-D MHD plasma-neutral gas simula- 
tions assuming invariance in the z-direction (d/dz ~ 0). 
In the simulations, a dense neutral gas cloud is moving in the 
y-direction impinging on a homogeneous magnetic field lying 
in the x-direction (see Fig. 

Here, we consider, for example, physical parameters for HVCs 
in the Magellanic Stream, i.e. the outer Galactic halo. For 
the high velocity clouds MSI- VI in the Magellanic Stream 
the particle density of the ambient plasma lies in the range 
of 10 



cm 



7 X 10 



-4cm-3 



(Weiner & Williams 1996 



Mirabel et al. |1979| ) for the Galactic halo and is even lower 

-6 



■at 



about 10~ cm~'^ - in the intergalactic medium. The temporal 
evolution of the plasma-neutral gas system is calculated by 
numerically integrating the following set of normalized fluid 
equat ions ( for a derivation of the fluid equations see Birk & 
Otto ( |1996| )): 

f--V.(p.) (2) 

(plasma continuity equation) 
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(3) 

( neutral gas continuity equation ) 



1, 



— {pv) = -V • (pvv) - -Vp +{V xB)xB 



dt 



d_ 

dt 

OB 

'dt 



-piyf^ {v - t>„) 



(PnVn) 



(4) 

(plasma momentum equation) 
1. 



-V • (pnVnVn) " ^^Pn - Pnl'21 (^n - 1^) (5) 

(neutral gas momentum equation) 

V x{vxB) (6) 

(induction equation) 
P 



dp ^ 



7pV • - 3 (7 - 1) ;/i2 ( P - — Pn ) (7) 

Pn 

(plasma pressure equation) 



dpn 

dt 



Vpn - TnPnV • Vn 
P 



-3 (7n - 1) 1^21 Pn 



(8) 



(neutral gas pressure equation) 
The quantities p, p, v, p^, p^, v^, B, e, 7, and 7,1 denote 
the plasma mass density, pressure, and velocity, the neutral 
gas mass density, pressure, and velocity, the magnetic field, 
the elementary charge, and the ratio of specific heats for 
the plasma and the neutral gas, respectively. The term vv 
represents the dyadic product {viVj)-j. The ratio of specific 
heats has been chosen as 7 = 7,1 = 5/3. 
In the presented simulations, we assume a continuous ioniza- 
tion equilibrium. So far, we are not interested in ionization 
phenomena, but in the dynamical evolution of the high velocity 
cloud, only. 

The frequencies 1/^2 ffi denote effective elastic collision 
frequencies. In order to guarantee momentum conservation, 
they have to fulfil the following relation (for details see Birk & 
Otto (|1996|)): 



^12P 



^tlPn 



(9) 

The symmetric collision frequency 1/^3 is a measure for the rel- 
ative collisional friction between the plasma and the neutral 
gas. Heuristically, we chose the reasonable ansatz 



'12 



^oPn 



rpl/2 



for the collision frequency i^?, where T 



(10) 

p/p is the normal- 
ized plasma temperature and j/q is a constant factor. 
All quantities X are normalized to typical values X of the sys- 
tem. Here we chose parameters typical for high velocity clouds 
in the Magellanic Stream and neutral gas clouds in the inter- 
galactic medium. Length scales are normalized to a typical ex- 
tension of a cloud of L = 100 pc 3 x 10^" cm. The particle 
density (ions and neutral atoms, respectively) is normalized to 
a value of n = 10^^ cm~'^. The magnetic field is normalized 
to ij = 3 pG. The actual magnetic field in the Galactic halo 
is not known. Upper limits from cosmic ray confinement, pul- 
sars, etc. indicat e field strengths of some microGauss (Vallee 



1990| , Beck et al. [1994j). We perform the simulation with an ini- 



findings show that even with this relatively weak initial field 
a substantial dynamical stabilization of the HVC results. We 
assume a monoatomic neutral gas, e.g. hydrogen atoms, such 
that the ion mass and the neutral atom mass are equal. Further- 
more, we assume only one species of single ionized ions. 
One may note that a pure hydrodynamic treatment is barely 
applicable to HVCs in the Galactic halo since the fluid approx- 
imation does not strictly hold. 

With the above choice for L, n, and B further normaliza- 
tions follow in a generic way. The mass densities p and pn 
are normalized to p = (mo + mi)n « min = mafi where 
nic, TOi, and rria, are the electron, ion, and neutral atom mass, 
respectively. Velocities are normalized to the Alfven veloc- 
ity VA = B/y/Airp, frequencies to the inverse Alfven time 
1/ta = WA /L, and pressures to p = B^/Stt. From the normal- 
ized equation of state for ideal gases p = nT the temperature 
of reference can be deduced as T = B^ /WnnkB where fee is 
the Boltzmann constant. 

For the above values of the system we find a mass density p w 
1.67 X 10^'^'' gcm^'^, the Alfven velocity va ~ 210 kms^^, 
the Alfven transit time ta ~ 1.31 x 10^^ s sa 416 000 yr, and 
the temperature T w 1.3 x 10'^ K. 

For numerical reasons, the density gradient of the modeled 
HVC is chosen to be smaller than a realistic one. Consequently, 
the neutral gas temperature in the neutral cloud is unrealisti- 
cally high. However, the dynamics are not significantly influ- 
enced by this choice since the mechanical forces are due rather 
to pressure than to density or temperature gradients. On the 
whole, we model the neutral gas as if its ionization energy was 
high enough to keep it neutral at lO^K. The ionization energy 
does not enter the calculation at any point, therefore being an 
arbitrary choice. 

The integration of the plasma-neutral gas equations is done 
on an equidistant 2-D grid by a second order leapfrog scheme 
where the partial derivatives are realized as finite differences by 
the FTCS method (Forward Time Centered Space). The conti- 
nuity, momentum and pressure equations are integrated with a 
flux correction of the source terms between the two half time 
steps while the induction equation is integrated by means of 
a semi-implicit algorithm (Dufort-Frankel ). Det ails of the nu- 
merical code can be found in Birk & Otto ( |1996| ). 
A schematic drawing of the simulation is given in Fig. ^. In or- 
der to simulate the infall of a compact dense neutral gas cloud 
perpendicular to a magnetic field in a hot plasma we first trans- 
form to the neutral gas rest frame. Then the initial configu- 
ration is defined by a dense neutral gas cloud centered at the 
origin {x, y) = (0, 0) of the system having the form 



Pn (r) 



(11) 



11111 I 1 

cosn ar 

Here, r gives the distance from the origin of the system while 
a = 1/2 defines the scale length for the neutral gas density 
gradient. The constants pnun and p^^ define the minimum 
and the maximum neutral gas density and are chosen to be 



Pn 



1 and 



25, respectively. The initial neutral 



tial field strength of one percent of the normalization field. Our 



gas density is shown in Fig. 
tially at rest, i.e., t;„ = 0. 



5|. The neutral gas cloud is ini- 
The neutral gas is assumed to 



6 



C. Konz, C. Briins & G.T. Birk: Dynamical Evolution of High Velocity Clouds in the Intergalactic Medium 




Figure 4. A high velocity cloud moving through a magnetized 
plasma 




Figures. The initial neutral gas density p„; figures are in non- 
dimensional units 



be initially in pressure equilibrium at a homogeneous pressure 
?5no = T^o (Pmin + Pno ) where Tnf, = 1 is the minimum tem- 
perature of the neutral gas cloud. The temperature profile is 
then given by Tn (r) — Pna/ Pn (r)(see Fig. ^). The neutral gas 
outside the cloud is about a factor 25 hotter than at the center 
of the cloud. This temperature gradient, however, is still small 
compared to realistic temperature gradients by factors of 10"^ 
- 10"* in typical HVCs. This is due to numerical reasons in so 
far as stronger gradients require an enhanced spatial and tem- 
poral resolution which makes a simulation extremely time- and 
memory-consuming. However, simulations with steeper gradi- 
ents in the temperature and density of the neutral gas compo- 




Figure 6. The initial temperature profile Tn for the neutral gas 

nent have shown qualitatively the same evolution of the system 
apart from the longer calculation time. Since we are interested 
in the longterm evolution of a plasma-neutral gas system we 
restrain ourselves to intermediate temperature and density gra- 
dients. 

The plasma itself initially has a homogeneous density pQ = 1 
and a similar temperature profile as the neutral gas: 

T(r)=ro-^^^. (12) 
cosh ar 

Here, Tq = 10 is the maximum temperature of the plasma. It is 
worthwhile mentioning that the external neutral gas is even hot- 
ter than the external plasma in this simulation. The neutral gas 
outside the cloud just serves as a reservoir of heat for the outer 
plasma which is heated up in the course of time via the colli- 
sional i^fj'term in the pressure equation. It is not dense enough 
to influence the motion of the plasma significantly due to the 
small value of the collision frequency = 8 x 10~^. Dur- 
ing the dynamical evolution, the neutral gas gathers the same 
velocity profile as the plasma motion. Initially, the plasma 
streams in the y-direction with the velocity profile 

^''^ 

where Vy^ ~ —0.1 is the minimum velocity of the plasma. 
This leaves the center of the neutral cloud at rest while plasma 
is streaming past the neutral cloud in the negative y-direction 
and across an initially homogeneous magnetic field in the 
x-direction: Bx = Bxo = 0.01, By — Bz — 0. Since the 
local Alfven velocity is given by ua — Bj y/Anp this plasma 
flow is superalfvenic at the beginning of the simulation due to 
the small magnetic field. This means that the plasma flow is 
able to compress and deform the magnetic field significantly. 
However, there is no shock to be found because of the high 
temperature Twl.SxlO^K and, connected with it, the high 
velocity of sound Cg = ^f^lrfk^Tjrai w 661 kms^^. 
The complete initial configuration is shown in Fig. ^. The 
neutral gas density is shown in logarithmic scaling while the 
lines indicate magnetic field lines and the arrows indicate the 
plasma flow. 
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Figure 7. The initial configuration for the neutral gas density, 
the plasma flow, and the magnetic field 



At the boundaries of the simulation box all quantities are 
extrapolated to the first order of the Taylor expansion. There 
are three exceptions, though. The plasma and neutral gas 
densities and the collision frequency 1^12 ^re chosen to be 
symmetric at both y-boundaries. Furthermore, the magnetic 
field component is chosen to be symmetric at the j/max- 
boundary because of numerical stability reasons. The most 
important exception, however, is given by the fact, that the 
plasma flow is maintained at the upper y-boundary by setting 
'Vy (j/max) — VyQ. Altogether, we get a plasma-neutral gas 
system where magnetic flux, plasma, neutral gas, and energy 
can freely cross the x-boundaries but where magnetic flux and 
mass flux are partly fixed at the j/-boundaries. 
The simulations presented here are all ideal in the sense that no 
resistivity apart from a small constant background resistivity 
rjo = 10^^ is applied due to numerical reasons. 
The calculations are done in a 2-D box with x going from —20 
to 20 and y going from —40 to 20. We use a uniform grid with 
103 grid points in the x- and 153 grid points in the j/-direction 
resulting in a constant grid spacing of 0.4. 
Since the system is treated ideally in a magnetohydrodynamic 
view the magnetic field lines are frozen into the plasma. 
Thus, they are deformed in the course of time by the plasma 
streaming past the neutral gas cloud. At the beginning of 
the simulation, both the neutral gas and the plasma at the 
origin are at rest. The collisional terms pi'f2 {"^ ~ ''^n) and 
PnJ^fi {vn — v) in the momentum equations lead to a deceler- 
ation of the plasma impinging upon the neutral gas cloud and 
to an acceleration of the neutral gas cloud itself. However, 
the acceleration of the neutral gas cloud is rather small due to 
the high neutral gas density of the cloud. Outside the cloud, 
plasma and neutral gas densities are comparable such that the 
neutral gas acquires the same velocity as the streaming plasma. 
The neutral gas flow is shown in Fig. ^. One finds a velocity 
pattern that resembles the initial velocity pattern apart from the 
fact that the center of the neutral gas cloud has moved towards 
the lower boundary. Furthermore, the plasma and neutral gas 



Figure 8. The neutral gas flow at i w 330ta 



is now flowing around the cloud as if it was a soft obstacle in 
a hydrodynamic flow. The center of the neutral gas cloud is 
now moving with a small velocity in the negative y-direction 
as a result of the continuous impact of plasma on the cloud. 
The exchange of momentum between the incoming plasma 
and the resting neutral gas cloud is performed via elastic 
ion-neutral gas collisions. The momentum exchange leads to 
an acceleration of the cloud in the initial rest frame, i.e., to 
a deceleration of the cloud in the plasma rest frame. This is 
consistent with ultraviolet absorption observations hinting at a 
deceleration of the HVCs as they approach the Galactic disk 
(Benjamin [T999[ Danly |l989| ). 
Fig. ^ shows the neutral gas velocity at the same time 




Figure 9. The neutral gas velocity at t « 330ta 

as the flow in Fig. || As already mentioned the center of the 
cloud is now moving with a small but finite velocity in the 
negative y-direction. 

The simulation time of about 330 Alfven times corresponds to 
a time interval of about 140 million years which is comparable 
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to crossing times of HVCs in the Magellanic Stream. 

A complete representation of the system after about 330 



Magnetic; Field 



magnetic barrier at the front of a HVC via the critical velocity 
effect (Konz et al. 2001 ). In turn, the formation of a magne- 




Figure 10. The neutral gas density, the plasma flow, and the 
magnetic field att ^ 330ta 



Alfven times is given in Fig. |10[ The arrows represent the 
flow of the plasma around the neutral gas cloud. Since the 
cloud does not act as a solid but as a soft obstacle there 
are no real stagnation points in front and behind the cloud. 
Nevertheless, the plasma flow encountering the cloud is partly 
redirected sideways similar to a hydrodynamic stream around 
a solid obstacle. Since the magnetic field lines are frozen in 
the plasma, magnetic flux is transported via the plasma flow 
towards the front edge of the cloud while the field lines on 
both sides of the cloud are draped around the cloud forming a 
magnetic barrier at the front and a magnetotail at the wake of 
the HVC. 

The formation of a magnetic barrier around the head of the 




Figure 11. The neutral gas density at t « 330rA 



neutral gas cloud resembles very much the formation of a 



totail is quite similar to the formation of a magnetotail by the 
well-studied interaction of the solar wind with unmagnetized 
Luhmann et al. 



planets (Luhman n [1986 
(Formisano et al 



1997) 



1982; McComas et al. 1987; Flammer et al. 



1992) and comets 



The neutral gas cloud itself is slightly deformed into a U-Uke 
structure, visible in Fig. |l^ in logarithmic scaling of the 
neutral gas density, together with a weak tail consisting of 
neutral gas being slowly stripped off the wings of the cloud. 
The stripping of neutral gas off the edges of the cloud results 
from the momentum transfer from the halo gas to the neutral 
gas particles. The stability of the cloud is not influenced 
significantly by this ongoing effect during the time period 
under consideration. Eventually, this abrasion of neutral gas 
forms a tail. The phenomenon of tail formation is consistent 
with the observation (see Fig. The upper column of 

this figure shows that the tail mainly consists of the warm 
component while the colder core of the neutral gas cloud 
remains unaffected by the stripping. After the formation of the 
magnetic barrier the resulting density profile of the cloud, as it 
is shown in Fig. [ll], remains almost constant in time as later 
time-steps of the simulation demonstrate. The HVC almost 
keeps its width and its maximum density. 
The dynamical stability of the neutral gas cloud is a direct ef- 
fect of the magnetic confinement by the accumulated magnetic 
barrier. This can be shown by performing the same simulation 
without any magnetic field. Figs. |l2| and [T^ show the neutral 
gas density at about 300rA and 540ta for this case. One finds 
the initially compact neutral gas cloud being first disrupted 
into two major parts and then completely falling apart. The 
cold neutral cloud is unstable against impact of hot plasma. 
The momentum and energy transfer from the incoming plasma 
to the front of the cloud leads to a disruption of the cloud into 
several parts. 

For comparison to the case of a magnetic field initially 
perpendicular to the flow direction. Figs. ^ and |l5] show 
the cloud and the magnetic field configuration for the case 
of an inclined magnetic field with initial field strengths 
Bxo = ByQ — 0.03. As in the previous simulation, the initial 
plasma inflow is purely in the negative y-direction. As long 
as there is a significantly strong component of the magnetic 
field perpendicular to the flow direction the formation of a 
magnetic barrier stabilizes the cloud in two ways. First of all, 
the increased magnetic field in the barrier exerts a magnetic 
pressure B'^/Stt on the cloud and partly compensates the ram 
pressure pv'^/2 of the plasma stream. Since the system is 
ideal in the magnetohydrodynamic description plasma cannot 
move across magnetic field lines. Therefore, the plasma flow 
is redirected by the magnetic barrier to stream around the 
HVC. On the inner side of the magnetic barrier, the magnetic 
pressure is compensated by the internal plasma pressure. 

As expected, we find a lower limit of the magnetic field 
strength sufficient to stabilize the HVC. For the parameters 
chosen, this limit is about 15nG meaning one order of 



C. Konz, C. Briins & G.T. Birk: Dynamical Evolution of High Velocity Clouds in the Intergalactic Medium 



9 




build-up of a magnetic barrier because the plasma pressure 
gradient would be even steeper due to the loss of plasma. 
That means that more magnetic flux is transported inside the 
magnetic barrier and the magnetic tail by the plasma heading 
for the cold center of the cloud. The magnetic pressure finally 
stops this plasma flow and isolates the plasma inside the cloud 
from the ambient halo plasma. In our simulation, the plasma 
inside the cloud takes up to 50% of the total density. This 
percentage, however, only reaches such a high level because of 
the relatively small density of the neutral gas cloud compared 
to the density of the ambient plasma. For a higher neutral gas 
density the plasma fraction inside the cloud only amounts to a 
few percent. 

Altogether, we have a system consisting of three different 



Figure 12. The neutral gas density at i w SOOta for the non- 




magnetic field 

Figure 13. The neutral gas density at i w 542ta for the non- 
magnetic case 



magnitude higher than the mean intergalactic magnetic field 
(Vallee 1983[ ). At the beginning of the simulation, the plasma 
density was set to be homogeneous. However, there is a 
pressure gradient due to the fact that the plasma inside the 
cloud region is much cooler than outside the cloud. Due to the 
finite temperature of a HVC and neutral gas ionization via UV 
photons, collisions and soft X-ray background there is always 
a small fraction of plasma inside a HVC amounting to the 
order of several percent of the total density. In the course of 
the simulation, this fraction is increased up to 50% by plasma 
streaming towards the center of the cloud along the field lines, 
i.e., along the tail and by a compression of the whole cloud due 
to the external plasma pressure. The plasma inside the cloud, 
however, is cooled by energy transfer to the cold neutral gas 
and would be subject to recombination events if recombination 
was included in the simulation. So far, recombination is not 
considered in our simulations. However, recombination of 
the cold plasma inside the cloud would further increase the 




Figure 15. The neutral gas density at t 
inflow magnetic field 



SOOta for an inclined 
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regions: an outer region of hot ambient plasma exerting the 
ram pressure on the region of the magnetic barrier and the 
inner region with cold plasma and neutral gas compensating 
the magnetic pressure from inside. Elastic collisions inside the 
cloud are much more frequent than outside since the higher 
neutral gas density overcompensates the lower temperature in 
Eq. (p^. Therefore, plasma and neutral gas inside the cloud 
have the same temperature and the same velocity. 
The magnetic barrier not only dynamically stabilizes the HVC 



k 




Figure 16. The magnetic field strength at f « SSOta 



but also thermally. The thermal conduction perpendicular to a 
magnetic field is drastically reduced compared to the thermal 



conduction parallel to the magnetic field (Braginskii 1965) 



The coefficient k± for the heat conduction perpendicular to 
the magnetic field decreases with increasing magn etic field 
strength as B^^ (Van der Linden & Goossens 1991 ). Fig. |l6| 
shows the magnetic field strength after about 330ta. At the 
front of the HVC, the magnetic field strength is increased 
by a factor 10, in the tail by a factor 4. This means that the 
perpendicular heat conduction is reduced by a factor 1/100 or 
1/16, respectively. In our simulation, heat conduction was not 
included so far Since the main part of the temperature gradient 
is perpendicular to the magnetic field lines no qualitative 
changes are to be expected if heat conduction is included. 

Heat conduction is suppressed and plasma cannot cross 
the magnetic barrier the temperature of the neutral gas cloud 
remains almost constant in time. Fig. |l7| shows the temperature 
profile of the neutral gas after about 330ta. Although a little 
bit deformed the initial temperature profile is almost main- 
tained. Steep temperature gradients are found in the region of 
the magnetic barrier surrounding the minimum temperature at 
the maximum of the neutral gas density. The temperature on 
the outside is still by a factor 10 higher than on the inside of 
the cloud. The magnetic barrier therefore insulates the cloud 
against heat exchange and heat conduction. 
In the tail, the temperature gradient is much smaller leading 
to a tail at an intermediate temperature. This relatively warm 




Figure 17. The neutral gas temperature at i sa 330ta 

neutral gas in the tail shows a higher velocity dispersion than 
the core of the cloud itself, a feature that has been confirmed 
by measuring the velocity dispersion both in the head and in 
the tail of HVCs (see Sect. 2). 
How is the neutral gas density in the simulation compared 



Keutrol Gas Density p„, tirae = 330.41 




Figure 18. A neutral gas density cut along the y-axis at t 
330ta 



to observations of head-tail structures? Fig. 18 shows a cut 
along the y-axis for the neutral gas density, i.e., along the axis 
of symmetry. In the plasma rest frame, the neutral gas cloud 
is moving to the right, i.e., in the positive y-direction. The 
density cut shows a steep gradient at the head of the cloud 
where the cloud is well confined by the magnetic barrier and 
a smooth almost exponential decay towards the tail. In this 
presentation, the HVC looks like a comet trailing in its wake a 
tail of material that is stripped off its surface in the course of 
time. For the HVC, the situation is quite similar Neutral gas is 
slowly stripped off the outer wings of the cloud by momentum 
transfer between plasma and neutral gas. The stripped gas 
atoms or molecules form a tail-like structure with a higher 
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temperature than the core of the cloud and a much smaller 
density. The tail itself is laterally shielded by magnetic field 
lines draped around the cloud forming a magneto-tail such that 
plasma and heat exchange between tail and ambient medium 
are reduced. The simulated HVC shows the same morphology 
and qualitatively similar density and temperature profiles as 
typical HVCs. Moreover, the velocity pattern (Fig. ^) shows 
the same characteristics as the observed HVC presented in 
Fig. ^: a fast moving cold core, slower moving warm boundary 
layer and a slowly moving tail. 

Finally, Fig. |l^ shows the column density for the neutral 
gas for a cut along the x-axis through the core of the cloud. 
Here, we calculated the column density by rotating the cloud 



i.e., a magnetic barrier forms similar to magnetopauses around 



Column density N^, 



330.41 




Figure 19. A neutral gas column density cut along the x-axis 
at t « 330rA 

around the y-axis and integrating along the line of sight lying 
in the z-direction. This corresponds to measuring the column 
density by integrating along the line of sight. Since we initially 
assumed invariance in the z-direction this is of course not exact 
but reasonable for clouds which are not too small. The result 
shows a smooth profile which is qualitatively very similar to 
column densities shown, for example, in Fig. |[ 
All in all, the simulation shows the importance of the magnetic 
field for the morphology of the HVC as well as for its stability. 



4. Discussion 

We presented the first plasma-neutral gas simulations of the 
interaction of HVCs with the Galactic halo gas. Previous nu- 



merical studies (Santillaan 1999 and Quilis 2001) focused on 
the dynamics in external. Galactic gravitational fields. The 
questions of thermal and magnetic insulation of the neutral gas 
clouds in the hot ambient medium as well as their stability were 
not addressed. Our two-fluid approach provides rather detailed 
information on the interaction region of the cold almost neu- 
tral HVC gas and the hot almost fully ionized ambient plasma. 
The Galactic magnetic field is locally draped around the HVC, 



comets (McComas 1987) and unmagnetized planets (Luhman 



1986). This barrier is of great importance for the thermal insu- 



lation of the cold, fast traveling cloud from the hot surround- 
ings. Furthermore, it prevents the plasma from diffusing into 
the neutral gas cloud. We conclude that the formation of mag- 
netic barriers contributes significantly to the dynamical stabil- 
ity of HVCs. Even for relatively weak initial magnetic fields in 
the Halo the magnetic barrier sufficiently stabilizes the cloud. 
In particular, unstable Kelvin-Helmholtz modes are not excited 
due to the stabilizing effect of a magnetic field component 
parallel to the flow. Thermal radiative losses that are not in- 
cluded in our contribution may further enhance the stability 
of the bo undaries against Kelvin-Helmholtz modes (Vietri et 
al. 1997 ). Within the regions of discontinuity and in the ad- 
jacent inner layers the physics becomes rather involved. In 
particular, some fraction of the stored magnetic field energy 
can be converted into heat by magnetic reconnection and thus 
can be responsible for the observed x-ray emission (Zimmer et 
al. 1997, Birk et al. 1998 ). Moreover, one may note that the ef- 
fect of ionization caused by the relative cloud halo movement 
(Konz et al. 2001) may enhance the formation of magnetic bar- 
riers by ionized boundary layers. 

A self-consistent three-dimensional simulation including all 
relevant physical effects seems to be highly desirable and is 
a promising task for the future. 

Finally, we note that the interaction scenario studied in this 
paper may also be of relevance in the context of cloudlets in the 
magnetospheres of Active Galactic Nuclei as well as in cluster 
cooling flows. 
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